Linear and non-Linear Instabilities in Unified Dark Energy Models 
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We revisit the paradigm of unified dark energy discussing in detail the averaging problem in this 
type of scenarios, highlighting the need for a full non-linear treatment. We also address the question 
of if and how models with one or several dark fluids can be observationally distinguished. Simpler 
and physically clearer derivations of some key results, most notably on the relation between the 
generalized Chaplygin gas and the standard (ACDM) 'concordance' model and on a Jeans-type 
small-scale instability of some coupled dark energy/dark matter models are presented. 



I. INTRODUCTION 

Observations have been steadily piling over the years 
that can only be 'explained', in the context of General 
Relativity, by the presence of so-called dark forms of en- 
ergy: ordinary baryonic matter and radiation are simply 
not enough [J 0, H, 0, [I| ■ In broad terms, a dark energy 
component (violating the strong energy condition) is re- 
quired to explain the recent acceleration of the Universe 
and a (cold) dark matter component to account for the 
amount of structure observed. 

Unfortunately, these components have never been ob- 
served directly. In fact, they may not exist at all. It 
is certainly possible to engineer dark energy and dark 
matter out of the picture, by modifying General Relativ- 
ity directly; often this process involves extra-dimensions 
with non-trivial topologies and/ or fiddling with TeVeS 
gravity sources @ (see also @,[K). At the moment, how- 
ever, most modifications seem far too ad hoc to be par- 
ticularly pleasing. In this paper, we will mostly work 
within the confines of General Relativity, i.e., we will as- 
sume that dark energy and dark matter are real entities. 
On the other hand, if this assumption is correct, it is 
truly mind boggling to find out that roughly 96% of the 
Universe should be in this dark form. Cosmology is thus 
hard-pressed to explain what these components are. 

Certainly the realization that dark energy and dark 
matter may not be independent entities should play a 
significant role. Although we frequently perceive them 
as being different, this does not have to be the case. The 
Chaplygin gas @, [K| , for instance, experienced a short 
burst of popularity precisely for being able to mimic both 
dark energy and dark matter (depending on the local den- 
sity), and still be just one form of energy — an exotic one, 
granted, but the same can be said of quintessence, k- 
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essence, and the rest of the phenomenological zoo. Thus, 
we should be open to the possibility that dark energy and 
dark matter are actually just different manifestations of 
the same underlying field, not different entities per se. 
On the theoretical front, the idea of a unified description 
of dark energy and dark matter has incredible heuristic 
potential and should be properly investigated. In case 
of failure, we would at least come back justified in treat- 
ing dark energy and dark matter as an independent pair. 
Either way, something of value could be gained. 

While the background and linear properties of unified 
dark energy models have been extensively studied by the 
community, the non-linear aspects of the problem have 
been largely ignored. In the following section, we will 
discuss in detail the averaging problem in the context of 
these models and show that non-linearities are a crucial 
aspect that cannot be neglected. 



II. AVERAGING IN UNIFIED DARK ENERGY 
MODELS 

The real Universe (henceforth denoted by manifold M.) 
displays hierarchical structures like stars, galaxies, clus- 
ters of galaxies, and so on; it is far from being smooth. 
The dynamics in M. is assumed to be entirely described 
by General Relativity (or some modified version of), as 
gravity is the relevant force at work on cosmological 
scales. In practice, however, due to the highly non-linear 
nature of the field equations, it is virtually impossible to 
solve them fully (even numerically) , except for a few high 
symmetry configurations; see, for instance, [Til ] . 

On the other hand, several observations indicate that 
the Universe looks increasingly smooth, as larger and 
larger scales are considered (typically over 100 Mpc). 
This average background (representing the global be- 
haviour of M) is routinely idealized as a completely fea- 
tureless manifold, hereafter denoted by manifold (M). 
The kinematics in (M) is essentially contained in the cos- 
mological principle, but what about its dynamics? This 
is actually a tricky question. Normally, it is assumed 
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that General Relativity applies just as well in (M) as it 
does in M. Yet, when we average M to obtain the back- 
ground, we are also averaging complex non-linear inter- 
actions. Averaging linear terms in the field equations is 
no big deal, but averaging non- linear terms is: this is be- 
cause they introduce backreaction terms and hence new 
dynamics. 

The common expectation however, is that this backre- 
action is negligible. Proving this, of course, is very hard; 
even checking it numerically is extremely difficult. As 
a matter of fact, 'averaging', in General Relativity, re- 
mains largely an unresolved problem [13 , [l3| . In connec- 
tion to this, it has even been suggested that the current 
acceleration of the Universe might be the result of these 
corrections, which would avoid the need for an exotic 
dark energy component (see for example [14J and refer- 
ences therein). As appealing as this may sound, however, 
such corrections are extremely hard to quantify, and cur- 
rently there is no convincing argument to show that they 
are sufficiently large to do the job. The few times that 
such an analysis has been attempted, by studying high 
symmetry configurations such as closely spaced sheets of 
matter separated by voids (which are arguably not a good 
approximation of the real universe), have shown them to 
be small, and hence incapable of producing an acceler- 
ated behaviour [ll| (see also [IB]). 

A related problem that will be relevant for what follows 
concerns the averaging of the energy content of the uni- 
verse described by the energy-momentum tensor of the 
matter fields. The evolution of (M) is usually assumed 
to be described by General Relativity with a source term 
given by the energy-momentum tensor of a perfect fluid 
with an energy density (e), where ( ) is an appropriate 
average, and a known equation of state, p = p(e). Even 
if we take for granted that for p — backreaction ef- 
fects are small, we can not ignore the fact that in general 
(p) p({ £ ))- An exception to this is the special case 
where the equation of state is given by p — we where w 
is a constant; then, we have (p) = w(e) = p((e)). How- 
ever, this is no longer true for the equation of state of 
the phenomcnological model known as the generalized 
Chaplygin gas \u[ . 
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where A and a are positive constants. Indeed, it is 
straightforward to check that 



(p) = -(A/e°')^-A/(sr=p((e)). 



(2) 



This is in fact a crucial problem in the analysis of all 
unified dark energy models not sufficiently close to a 
ACDM model (which for t he g eneralized Chaplygin gas 
is characterized by a = [lj|]) for the above effects to 
be negligible indeed, this invalidates quite a few sim- 
ple analyses. In fact, it has been shown that in some 
models the transition from matter to dark energy domi- 
nation may never occur or may do so much earlier than 



expected in the absence of perturbations Of 
course, this problem only occurs when perturbations be- 
come non-linear on certain scales. If, on the other hand, 
perturbations are linear on every scale then 



(p) = <p(<e>(l + *))>~<P(<e>)>, 
if 5= (e- (e))/(e) < 1. 



(3) 



III. SCALAR FIELDS 



In this paper we will use a (— 
Now, let us consider the action 



-) metric signature. 



d x ^—gC(X, i 



(4) 



where C is the Lagrangian density, is a real scalar field 
and 



(5) 



is the kinetic term. As a simple example, consider the 
case of a classical scalar field <j> governed by the so-called 
canonical Lagrangian 



c = x - v{4>) , 



(6) 



where V(<f>) is some scalar potential. The energy- 
momentum tensor may be obtained by varying the action 
relative to the inverse metric, which yields 

T^{4>) = V^V„0 + (X - V(<f>)) gitu . (7) 

If we now make the following identifications 

u„ = ^dL, e = X + V{(j>), p = X-V(<j>), (8) 



and plug them into (|7|) we obtain 

= (e + p)u^u v + pg» v ; 



(9) 



hence the energy-momentum tensor of <j> can be written 
in the form of a perfect fluid. Recall that in ©, u M is the 
4-velocity field describing the motion of the fluid (if V M </> 
is timelike) while e and p, are its proper energy density 
and pressure, respectively. 

The above is easily generalized for an arbitrary La- 
grangian density of the form C(X,cf>) [2l[. This defines 
a new class of scalars broad enough for our purposes. 
Again, we can still explicitly rewrite the energy momen- 
tum tensor for this model in a perfect fluid form, by 
means of the following identifications 



e = 2Xp. x -p, p = L(X,<f>). (10) 



From this it follows that if p = p(X), then e = e(X). 
Unfortunately, it is not always possible to invert s(X), 



3 



and obtain X{e) but when it is, the fluid has an explicit 
isentropic equation of state p — p{s). A useful example is 
p oc X n which, as one can easily check, describes a con- 
stant u> — l/(2n — 1) fluid. In particular, when n = the 
scalar corresponds to a standard cosmological constant, 
n = 1, to a massless scalar field, n = 2, to background ra- 
diation, and so on. In the limit of large n, the scalar field 
can be interpreted as dust (a pressureless non-relativistic 
fluid). 

How about the generalized Chaplygin gas? We can 
retro-engineer IjlOp to ascertain the necessary Lagrangian 
by solving 



2Xa 



A e 



x 



A 



(11) 



a non-linear differential equation. Dividing by e and 
rewriting everything in terms of a new function £ = 
A/e 1+a , we end up with a much nicer linear version 



1 = -- 



l + a 

which is easy to solve. The solution is 
£ = 1 - (2X)^ 



(12) 



(13) 



and thus the Lagrangian which reproduces the general- 
ized Chaplygin becomes 

C{X)=p(X) = ~ = -A^^ , (14) 

with < 2X < 1. 

Finally, the sound speed for any scalar field (canonic or 
otherwise) is also well defined and determined in linear 
theory to be [21 1 



2 _ P,x 



s,x £,x + 2XC t xx 



(15) 



contrast this to the isentropic result Sp/5e. It follows that 
quintessence has a constant sound speed of unity. This 
is another reason why quintessence cannot reproduce the 
generalized Chaplygin gas. On the other hand, if we 
apply this formula to the scalar field governed by (fT?)) . 
the sound speed does come out the same as the isentropic 
one and all is well. 



IV. INTERLUDE: THE ACDM LIMIT 

Implementing the generalized Chaplygin gas as a scalar 
field obeying a particular Lagrangian is useful on a num- 
ber of levels. For one, it allows us to prove that the 
a = limit of the generalized Chaplygin gas is totally 
equivalent to an ordinary ACDM model. This fact plays 
an important role when comparing unified dark energy 
models with observations. While obviously true in the 
absence of perturbations, the need to explicitly demon- 
strate this beyond a zero-order equivalence, only became 



apparent after [22j appeared. In it the surprising claim 
was made that the linear evolution of perturbations for 
ACDM and the a = generalized Chaplygin gas model 
actually differed. Nevertheless, in it was shown that 
this was not true (see also (23|; their equivalence to first 
order in the metric perturbations was established and it 
was equally argued that the correspondence went well 
beyond linear order. Here we give a much simpler proof 
of this equivalence. If we expand (fT4|) in a power series 
and take the limiting case a — > we obtain 



p(X) = lim -A 1/1+a 
= -A + Q- {2X) C 



1 - 



a , . i+o 
(2X)^ 



l + a 
-A. 



(16) 



Thus, everywhere in M.. the Lagrangian decomposes 
nicely into a cosmological constant plus 'matter', thus 
demonstrating the equivalence between a = and 
ACDM to any order; gravity alone does not distinguish 
between the two. Moreover, given that both are, to a 
certain extent, simply toy- models with somewhat shaky 
motivations from fundamental physics, this is probably 
as far as they can be meaningfully compared. 



V. A TRULY 'ATOMIC FLUID? 

We have shown that the a = limit of the general- 
ized Chaplygin gas is equivalent to ACDM. An obvious 
follow-up question is whether this equivalence between a 
single unified dark energy fluid and some family of min- 
imally coupled components is valid in general. Recently 
[HI has argued that it is always possible to split a single 
unified dark energy fluid into different minimally coupled 
components or to combine several fluids into a single fluid 
behaving in exactly the same way as the original mixture, 
from a cosmological point of view. Although this is ob- 
vious in the absence of perturbations, it was also argued 
that this degeneracy went beyond the background level. 

Before trying to answer the above question, we need to 
clearly define what we mean by a single fluid. It is clear 
that by allowing the complexity of the fluid to be arbi- 
trarily large, for example by considering very high order 
tensor fields, we may in principle get disproportionately 
non-trivial dynamics. However, in this case, do we have 
one fluid or several minimally coupled fluids? In gen- 
eral, a complex single fluid description is also prone to 
a multi-fluid interpretation. In this paper, we shall refer 
to a single indivisible ('atomic') fluid as one whose La- 
grangian density is a function of a real scalar field which 
has the form C(X,<j>). We shall focus our attention on 
the isentropic subfamily characterized by an equation of 
state of the form p = p(s). 

It goes without saying that the energy-momentum ten- 
sor of any fluid can always be split into several pieces. In 
fact there is an infinite number of ways to accomplish 
this. The critical question, however, is how to interpret 
any such decomposition. Are the resulting parts real 
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physical fluids? Do they exist independently from one 
other? We have already seen that the answer is no if the 
generalized Chaplygin gas does originate from the single 
scalar field governed by (Tl4"|) . Each piece is then a virtual 
component, i.e., without independent existence (except 
in the special a = case). In fact, the evolution of in- 
dividual virtual parts is not, in general, constrained by 
causality. Conversely, if we are in the presence of various 
fluids we can also add their energy- momentum tensors. 
However, the dynamics of the resulting fluid could be 
very complex and in general it will not be describable by 
a real scalar field with a Lagrangian density of the form 
L{X,4>). 

From a cosmological point of view all the relevant in- 
formation is contained in the energy momentum tensor 
which acts as a source for the gravitational field. Pend- 
ing laboratory evidence (which in principle can detect not 
only fields themselves but even couplings between them) 
we are only sensitive to the total energy-momentum ten- 
sor and consequently [24| argues that cosmology alone 
does not provide useful information on whether a single 
unified dark energy fluid or a family of minimally coupled 
fluids is responsible for the observations. 

However, if we have a single fluid described by an ar- 
bitrary equation of state p = p(e) and consider the evo- 
lution of very large wavelength perturbations in a homo- 
geneous and isotropic universe we know that, given local 
initial conditions for H and e, the subsequent evolution 
is the same in any such patch. This is no longer true if 
we consider a family of minimally coupled fluids. Let us 
consider large wavelength flat patches characterized by 
the same value of H (flatness implies that the value of e 
is also the same in any such patches). It is then easy to 
show that we may find any number of families of three 
minimally coupled fluids with the same initial conditions 
for e, e, p, and p but with very different subsequent evo- 
lutions. As an example consider the two families of three 
fluids characterized by a constant w with 

• Wi = — 1 = —W3 and W2 = 0, with £j/e = 1/3 for 
i = l,2,3 

• wi = — V6/3 = — W3, with Sj/e = 1/2 for j = 1,3 
and £2/^ = 

We see that if we are in the presence of various fluids, a 
similar behaviour at a given time may well lead to very 
different behaviours at a later time. Moreover, we are 
even allowed to decide how much freedom we wish to 
give ourselves: we may set up adiabatic or iso-curvature 
fluctuations (or indeed a combination of them). Clearly, 
the former situation will be more restrictive than the lat- 
ter. 



VI. LINEAR INSTABILITIES IN UNIFIED 
DARK ENERGY MODELS 



sound speed c 2 = but this is no longer the case in the 
context of more general models. It is easy to show (2f| 
that the evolution of the sound speed in the absence of 
non-linear perturbations is fully determined by the back- 
ground evolution and is equal to 



dp 
dp 



1 d 
3HdH 



H 



1 

q -2 



(17) 



where q = — a/(aH 2 ) is the usual deceleration parameter. 
The sign of c 2 s is thus linked to the background dynamics 
and in particular to how fast the transition from dark 
matter to dark energy occurs. 

A simple illustration of the above point is provided by 
the generalized Chaplygin gas equation of state given by 
Eqn. (JTJ, which corresponds to 



-A/e 



l+a . 



(18) 



phenomenologically we allow a to take any value. For 
a > this is the generalized Chaplygin gas itself, for 
a = this is equivalent to ACDM and for a = — 1 and 
A < 1 it is equivalent (in the absence of perturbations 
only) to a standard quintessence model with constant w. 
We shall not consider the case with a < — 1 since this 
would imply a universe which is dark energy dominated 
at early times and we are only interested in the late time 
behaviour of dark energy. This fluid illustrates the above 
points. For a > —1 its energy density evolves as 



£0 



n n \ 3(l+a) 

v a 



1/1+a 



(19) 



where A = A/eQ +a is a constant. The sound speed has 
the form c 2 s — —aw and consequently c 2 s > for a > 0. 
We clearly see that the sign is related to the speed of the 
transition from the matter dominated era to the dark 
energy dominated era. If it is steep enough (specifically, 
faster than in ACDM), c 2 will be positive, otherwise it 
will be negative. 

On the other hand, consider a universe with matter 
and a standard scalar field with the Lagrangian of Eqn. 
(|6|) and consider what happens if we treat cold dark mat- 
ter (CDM) and dark energy (DE) as a single isentropic 
fluid. This may be justifiable, to a certain extent, if there 
is a strong coupling between the dark matter and dark 
energy sectors as is the case in (2|| H3, HH ■ In the absence 
of perturbations, the energy density of this effective fluid 
satisfies 



- = -3H(l+u), 

£ 



(20) 



with w = pde/(£de + £cdm) > — 1- Since e > we 
obviously have i < 0, and the sign of c 2 only depends on 
the pressure evolution which is given by 



The simplest unified dark energy model (the general- 
ized Chaplygin gas in the a = limit) has a vanishing 
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where the slow-roll approximation has been used. Hence, 
in this case 

c? = f<0. (22) 

This instability [25[ will occur in any case where p > 
(assuming i < 0). A particular example of this is the 
Jeans-type instability studied in [H, 113, H, H, HJ, 
which occurs for models where there is a sufficiently 
strong coupling between dark matter and dark energy. 
(Note, however, that such large couplings are strongly 
constrained by several equivalence principle type experi- 
ments, so the cosmological relevance of these models re- 
mains unclear.) If wde is a constant, it is straightforward 
to show that 

c 2 _ w DE (l + Wde)£de ^3) 
(1 + wde)£de + £CDM 

Interestingly, the models where this instability is absent 
are precisely those where one can get wde < — 1 (as long 
as w DE > -1 - ecdm/sde) HE m M, H, HI- Hence, 
the above argument for c 2 s < no longer holds if one 
allows wde < — 1. 

VII. CONCLUSIONS 

In the so-called concordance cosmology model, a range 
of observational data is used to postulate the existence 



of two dark fluids for which so far there is no direct ex- 
perimental evidence. It is therefore natural to consider 
scenarios where these fluids are coupled, or even different 
manifestations of a single fluid. 

In this paper, we have re-visited unified dark energy 
models and discussed some of their distinguishing fea- 
tures. In particular, we have addressed in some detail 
the averaging problem, showing that non-linear effects 
cannot, in general, be ignored. We have also clarified the 
relation between a particular class of UDE models and 
the standard ACDM paradigm, and the physical reason 
behind a Jeans-type instability that has been identified 
in sufficiently coupled models. 

We have equally discussed the question of whether 
models with one or several dark fluids can be observation- 
ally distinguished showing that in the latter case, similar 
initial conditions for the total energy density and pres- 
sure can lead to very different outcomes. This degeneracy 
is absent in single fluid models. 

Finally, we would like to point out that a laboratory 
detection of one or more dark components is not the only 
way of ultimately breaking this degeneracy. Equivalence 
principle tests, for example, can provide key information 
on the dark sector. Even more promising is the prospect 
of using astrophysical measurements of varying couplings 
to reconstruct both the dark energy equation of state and 
a measure of the dark sector's coupling to the standard 
model as a function of redshift [3l| . This will be further 
explored in future work. 
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